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The zero temperature response of an interacting electron liquid to a time-dependent vector potential 
of wave vector q and frequency u such that q <C qp, qvF Cw< Ep/h (where qp, vp and Ep are 
the Fermi wavevector, velocity and energy respectively) is equivalent to that of a continuous elastic 
medium with nonvanishing shear modulus fi, bulk modulus K, and viscosity coefficients n and f . We 
establish the relationship between the visco-elastic coefficients and the long-wavelength limit of the 
"dynamical local field factors" Gi(T)(g, id), which are widely used to describe exchange-correlation 
effects in electron liquids. We present several exact results for (i, including its expression in terms of 
Landau parameters, and practical approximate formulas for /i, 77 and £ as functions of density. These 
are used to discuss the possibility of a transverse collective mode in the electron liquid at sufficiently 
low density. Finally, we consider impurity scattering and/or quasiparticle collisions at non zero 
temperature. Treating these effects in the relaxation time (r) approximation, explicit expressions 
are derived for /_t and r\ as functions of frequency. These formulas exhibit a crossover from the 
collisional regime (lit <C f), where fi ~ and n ~ nEpr, to the collisionless regime {lot 3> 1), where 
/i ~ nEp and r\ ~ 0. 
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I. INTRODUCTION 

The response of a solid body to external macroscopic 
forces is described by the theory of elasticity.tl In a ho- 
mogeneous and isotropic bodyo the response is controlled 
by two real elastic constants, the bulk modulus K and 
the shear modulus fx; dissipation is negligible. 

The macroscopic response of a liquid system, on the 
other hand, is usually described in terms of the Navier- 
Stokes equationlj of classical hydrodynamics. This is, at 
first sight, very different from elasticity. First of all, by 
the very definition of liquid, the shear modulus vanishes. 
Second, there is dissipation, due to the two viscosity co- 
efficients r\ and £ - the "shear" and "bulk" viscosities 
respectively. Only the bulk modulus remains approxi- 
mately the same in the liquid as in the solid state. 

Such a sharp distinction disappears at finite fre- 
quencies, where liquids develop solid-like characteristic, 
namely, a non vanishing shear modulus. Both liquids and 
solids follow a common visco-elastic behavior, which can 
be mathematically described by a single set of equations 
(say the equations of elasticity) with complex frequency- 
dependent elastic constants 



and 



K(lo) = K(lu) - iwC(w) 



fb(ui) — h(uj) — iujrj(w) . 



(1.1) 



(1.2) 



The visco-elastic coefficients K, /j,, £, rj on the right hand 
side are all real functions of frequency. 



The crucial parameter that controls the prevalence of 
solid-like or liquid-like behavior in the liquid is lot, the 
ratio of the frequency to the inverse of the relaxation 
time t - the time it takes the system to return to ther- 
mal equilibrium after being slightly disturbed from it. 
If lot -C 1 one is in the collision-dominated (or hydrody- 
namic) regime, in which /j.(lo) is negligible and rj(u)), C( w ) 
are finite. If, on the other hand, lot 3> 1 one is in the 
collisionless (or elastic) regime, where fi(u>) has a finite 
value, while the viscosities are small. In either case, the 
bulk modulus does not show a significant dependence on 
frequency. 

In this paper we explore the possibility of describing; 
the long wavelength dynamics of a quantum Fermi liquidcl 
near the absolute zero of temperature in terms of classical 
visco-elastic equations of motion. Limiting ourselves to 
the linear response of the quantum liquid to an external 
vector potential A(q, u>) of wave vector q and frequency 10 
we shall show that the visco-elastic description is possible 
(and useful) in the regime 



q<.qp, q < oj/vf, 



(1.3) 



where vp is the quasiparticle Fermi velocity and qF is 
the Fermi wave vector. In other words, the frequency 
must be high compared to the characteristic energy of 
quasiparticle-quasihole pairs at wave vector q, which 
tends to zero when q — > (see Figure [j]). 

The viscoelastic coefficients will be expressed in terms 
of the long wavelength limit of the dynamical local 
field factors G/,(t)(?|W): these are mathematical con- 
structs (defined below) that are widely used to describe 
exchange-correlation effects in Fermi liquids. 
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cj = qv F 



FIG. 1. The region below the line to — vfQ is the quasi- 
particle-quasihole regime. The viscoelastic approach applies 
to the region uj 3> Vpq. The line u> = 1/t separates the col- 
lisional viscoelastic regime and the collisionless viscoelastic 
one. 



Most of this paper is devoted to the task of calculating 
the viscoelastic coefficients of an electron liquid (both in 
three and two dimensions) in the limit of to — > 0, that 
is, in practice, for ui <C Ep but still satisfying the con- 
dition (1.3). Such coefficients are particularly relevant 
in the framework of Time-Dependent Density Functional 
Theoryjj where they fully determine the low-frequency 
regime. 

We first consider the case of a uniform (translationally 
invariant) electron gas at the absolute zero of temper- 
ature. In this case the low-frequency elastic constants 
K and \i can be expressed exactly in terms of the Lan- 
dau parameters Fq, F\ and F% - at least insofar as the 
Landau theory of Fermi liquids applies. The result for 



the bulk modulus [K = n d ' e(n)/dn , where e(n) is the 
ground-state energy density anxljji is the particle density] 
has been known for a long timdSQ and can be straightfor- 
wardly ressaluated from the knowledge of the ground state 
energy.Qfl The result for the shear modulus is (to the best 
of our knowledge) new and, unfortunately, not so easy to 
evaluate. For this reason, we propose an approach "a la 
Wigner" , namely we calculate the shear modulus at both 
high and low densities - where the calculation can be 
done with relative ease - and interpolate between these 
two limits. The proposed interpolation function is close 
to the results of recent mode-coupling calculations of the 
dynamical local field factora at sufficiently high density. 

We proceed in a similar way to the calculation of the 
viscosities. First the shear viscosity r] is analytically cal- 
culated in the high density limit, making use the for- 
malism of Nifosf et al.p for the imaginary part of the 



dynamical local field factor. Then we devise a numerical 
fit that reduces to the analytical result in the high den- 
sity, limit and reproduces the numerical data of Nifosf et 
al.p at lower density. The bulk viscosity £ is found to be 
approximately zero in this approach. 

The last part of the paper is devoted to the treatment 
of relaxation effects caused either by collisions with im- 
purities or by collisions between thermally excited quasi- 
particles. We assume that both effects can be described 
by a single relaxation time t, such that 1/t ^ Ep, and 
make use ofi-the standard "relaxation time approxima- 
tion" (RTA)ll2l to approximate the collision integral in 
the quasiparticle transport equation. The low-frequency 
regime now splits into two distinct regimes: collisional 
(lot <C 1) and collisionless (lot 3> 1). The restriction 
given by Eq. (1.3) remains in force in both regimes. By 
solving the transport equation in the RTA we obtain ex- 
plicit expressions for the elastic and viscous coefficients. 
In the collisionless regime, these expressions reduce to the 
ones derived in the first part of the paper. In the colli- 
sional regime, they are very different: the shear modulus 
vanishes (as expected for an ordinary liquid), and the 
shear viscosity tends to the limit rj = /it where /x is the 
collisionless shear modulus. Our simple analytic expres- 
sions clearly exhibit the-crossover from the collisional to 
the collisionless regime.EJ 

W e have emphasized the importance of the condition 
(1.3) that assures the possibility of a visco-elastic descrip- 
tion of the dynamics of the Fermi liquid. What happens 
if this condition is violated? The behavior of the micro- 
scopic current-current response function of a Fermi liq- 
uid changes dramatically as one goes from the q <C io/vp 
regime to the q 3> uj/vp regime, even though q and ui 
remain small compared to qp and Ep respectively. The 
physical reason is that the response in the second region 
is dominated by electron-hole excitations that are absent 
in the first. Because of the change in the character of the 
response, the current does not obey classical visco-elastic 
equations of motion in the second regime. Alternatively, 
if one insisted on casting the equation for the current in a 
visco-elastic form, one would be forced to use visco-elastic 
coefficients that diverge in the q — ► limit. This shows 
that the visco-elastic theory is not a natural description 
of the physics for to < qvp. 

This paper is organized as follows: 

In Section [ij we briefly review elasticity, hydrody- 
namics, and the local field factor representation of the 
current-current response functions of a Fermi liquid. We 
establish the relationship between the dynamical local 
field factors and t he frequ ency-dependent visco-elastic 
coefficients of Eos. (QQ. 

In Section [II we derive an exact expression for the 
shear modulus of the Fermi liquid at T = in term of 
Landau parameters, and a rigorous upper bound on the 
value of the elastic constants. 



In Section [V we present approximate analytical ex- 
pressions for the evaluation of the visco-elastic coeffi- 
cients of an interacting electron liquid as functions of 
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density. These expressions are used to discuss the pos- 
sibility of a transverse sound mode in the low density 
electron gas. 

In Section ^ we include electron-impurity and ther- 
mally induced quasiparticle collisions via Mermin's relax- 
ation time approximation. We provide explicit formulas 
for the frequency dependent (on the scale of the inverse 
relaxation time) shear modulus and viscosity exhibiting 
the crossover between the collisional and the collisionless 
regimes. 



II. VISCOELASTIC CONSTANTS OF A FERMI 
LIQUID 



Both the current and the vector potential can be written 
as sums of longitudinal and transverse components (par- 
allel and perpendicular to q respectively) j = jh + Jt, 
A = Al + At and the equations of motion for longitu- 
dinal a nd t ransverse components decouple. The solution 
of Eq. (U) is 



n/m 



+2(1-1) ^] 



mu)- 



-A L (q,co) (2.7) 



and 



3T{q,0J) = z „ —i A T (q,u). 



~n? mu' 2 



(2.8) 



A. Macroscopic equations 

The equation of motion for the elastic displacomcnt 
field u(r, t) in a homogeneous and isotropic solid isEJ 



d 2 u 



K 



V(V • u) + fiV 2 u + F(f,t) 
(2.1) 



where K and /x are constants, known as the bulk and 
the shear modulus respectively, d is the number of space 
dimensions, F{r,t) is an externally applied volume force 
density, n is the equilibrium number density, and m is 
the mass of the particles. 

We consider periodic forces of the form 



F(r,t) = ^ujef^ + tc, 
which induce periodic displacements 

u(r,t) = u{q,u)e i{ -* f -^ + c.c. 



(2.2) 



(2.3) 



In order to make contact, later, with microscopic theo- 
ries of Fermi liquids, we wcite the force as the time deriva- 
tive of a vector potentialO 



F(f, t) = n 



dA(r, t) 
dt 



and introduce the current density 

du (r, t) 



dt 



(2.4) 



(2.5) 



as its conjugate field. The equation of motion (2.1), writ- 
ten in terms of the Fourier transform of the current den- 
sity, takes the form 



~K 




— 


n 




n 



q[q-j(q,u)] 



--— j(q,v) - iumA(q, u). (2.6) 
n iuj 



We note that ji(g, lu) is related to the induced density 
change ni(q, uj) by the continuity equation 



ni(<f,w) 



7l0?>) 



(2.9) 



and the longitudinal vector potential Ai,(q, to) is equiva- 
lent, modulo a gauge transformation, to a scalar potential 
V(q, u>) such that 



A L (q,Lo) = ^V(q,u). 
u> 



(2.10) 



In writing these equations we have assumed that the 
external field Al includes, self-consistently, the contri- 
bution of the mean electrostatic field generated by the 
density fluctuation n\ [the Hartrce field Vh — v(q)ni\. 

Let us now consider the classical hydrodynamical equa- 
tion for the current density in a liquid.0 In the linear 
approximation with respect to j and A it has the form 



iojmj(q,uj) 



q[q-j(q,ui)] 



luirj q -> _ -> 

— - iumA(q,u), (2.11) 

n ioj 



where we have used the continuity equation ( |2.9| ) 
to rewrite the hydrostatic pressure term in the Euler 
equationQ as 



- Vp(g, lo) 



dp(n) _± dp(n) 
-\/ni(q,(J) = — q 



dn 



dn 



q-j(q,u) 



IU! 



(2.12) 



The constants rj and £ are the shear and bulk viscosity 
coefficients respectively; p(n) is the equilibrium pressure 
as a function of density and is related to the bulk modulus 
by the relation K — ndp(n)/dn. 

Equations (2.11) and ( |2.6| ) are very similar. The 
hydrodynamic equation differs from the elastic equation 
through the following replacements: (i) The shear modu- 
lus jj, is replaced by the imaginary quantity —iLurj, which 
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vanishes at uj — in agreement with the notion that a 
liquid has no resistance to shear, (ii) The bulk modulus 
K acquires an imaginary part —iuj^ where £ is the bulk 
viscosity. 

These observations suggest the use of a single language 
- say that of elasticity theory - to describe both the liquid 
and the solid. In this generalized scheme, the equation 
of motion becomes 



iumj(q,uj) = 



AM 



9l<f- J(<fM] 



jliui) q 2 -> 

—j{q,u) - wnA(q,u), (2.13) 



where fl(uj) and K(uj) are the frequency-dependent vi sco- 
clastic constants defined in the introduction [Eqs. 

O)]- 

The solution of this equation of motion is 



j L (q,uj) = 



n/m 



1 - 



n z \ a) n 



rnuj 2 



-A L (q,io) 



and 



jT(<?M 



n/m 



-A T (q,uj) 



(2.14) 



(2.15) 



The key difference between a solid and a liquid is that 
the solid has an essentially real fl (u finite, n ~ 0), 
whereas a liquid has an essentially imaginary fi (rj finite, 
[i ~ 0). As for the generalized bulk modulus, its real 
part (related to the compressibility) is nearly the same 
in the two phases. The bulk viscosity is generally of the 
same order of magnitude as the shear viscosityo but, in 
the case of the electron liquid, it will be shown to vanish 
within the mode-coupling approximation of Ref. |9| 



B. Connection with microscopic linear response 
theory 

Let us now turn to the microscopic formulation of the 
linear response of a homogeneous, isotropic body, sub- 
jected to an external vector potential A(q, uj). The proper 
response functions xl(<1,u) and xr(<f, (longitudinal 
and transverse respectively) are defined by the relations 

jL(T)(q,w) = XL(T){q,oj)A L{T) (q,uj). (2.16) 

A useful way of representing Xl(t) 

X°(T)(g» 

1 + v{q)G L(T) (q, uj)(q 2 /u 2 )x° L{T) (q, uj) ' 

(2.17) 



XL(T)(q,u) 



where 



n 
rn 



E 

k 



f(e k 



k+q] 



W + tk+q — £fe 



(2.18) 



are the longitudinal (transverse) response functions of 
the noninteracting electron gas, = k 2 /2m is the free 
particle energy, f(ek) is the Fermi distribution function, 
and fci,(T) is the longitudinal (transverse) component of 

k relative to q, and v (q) is the Fourier transform of the 
interaction [v(q) = Aire 2 /q 2 and v(q) = 2iTe 2 /q in three 
and two dimensions respectively]. The dynamical local 
field factors G^^iq,^) are effectively defined by these 
equations: they take into account exchange-correlation 
effects beyond the random phase approximation. 

Note that, with these definitions, the longitudinal local 
field factor Gl coincides with the more familiar scalar 
local field factor G used in the theory of the density- 
density response function in Ref. |l^ for example. 

Let us now consider the long wavelength limit (qvp <C 
uj) of Xl(t) (<Zj w )- Expanding the noninteracting response 
functions (2.18) we find 



XL(T) (<7>) 



1 



(d) E F nq 2 



l L{T) 



cL(T) 



M 



no- 



where 



fxcL(T){u) = - ]imv(q)G L r T )(q,uj) 



(2.19) 



(2.20) 



are complex functions of frequency, which satisfy 
Kramers-Kronig dispersion relations between their real 
and imaginary parts (see discussion in the next section) 



and = 6/5, a 



3/2, o4 3) = 2/5, a 



(2) 



1/2. 



Compari ng Eq . (2.19 ) wit h the macroscopic response 
functions ( |2.14| ) and (|2.15|) we are led to the following 
identifications: 



AM 



a^nE F + n 2 f xcT (u>) 



(2.21) 



and 



K(u>) 



a 



(d) 



„(<*) 



nE F 

/xctM 



(2.22) 



Separating the real and imaginary parts of these equa- 
tions, and taking the limit uj — > (but still with uj 3> 
qvp), we arrive at the promised expressions for the elastic 
and viscosity coefficients in terms of the long wavelength 
limit of the local field factors: 



offlnEp 



n 2 Re/ xcT (0), 



(2.23) 
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K 



(<*) 



a 



-n 2 Re 



2-- d )a T 



uEf 



/xci(O) - ( 2 - - ) / xcT (0) 



(2.24) 



?? = — n lim 



C = -n 2 lim 



/Im/ xcL (o 



o V 



2\ lmf xcT (uj) 
d I ui 



(2.26) 



Equations ( ^.23|^.26[ ) are the main result of this Sec- 
tion. We underline the fact that they have been obtained 
under the assumption qvF <C oj. Outside this regime, 
e.g., for lu < qvp, the microscopic response functions 
do not yield visco-elastic equations of motion or, equiv- 
alently, the visco-elastic coefficients diverge for q — ► 0. 
The analysis of this non-viscoelastic regime is beyond 
the scope of this paper. 



III. THE CALCULATION OF THE ELASTIC 
CONSTANTS: RIGOROUS RESULTS 

A. Expression in terms of Landau parameters 

The elastic constants of a Fermi liquid can be exactly 
expressed in terms of Landau parameters, insofar as the 
Landau theory of Fermi liquids is valid. To see this, we 
begin by deriving the equation of motion for the quasi- 
particle distribution function in the presence of a slowly 
varying vector potential A{r 1 1). n 

Following the discussion of Nozieres and Pinesta we 
treat the quasiparticles as a gas of non-interacting clas- 
sical particles governed by a self-consistent Hamiltonian 

H qP (r,p) = ^+A(r,t) + fpp' Sn p'(r, t) 
p' 

p-A{r,t) 



<P) 



nr 



^2fp P >5n p ,(r,t), (3.1) 



where ep is the quasiparticle energy, V ' p ep : 



p/m* 



is the quasiparticle velocity, to* is the effective mass, and 
f pp > are Landau parameters. This approach is justified 
in the limit of zero temperature and zero frequency, since 
the excited quasiparticles have an essentially infinite life- 
time in this regime, and their mutual collisions are neg- 
ligible. 

The self-consistent nature of the quasiparticle Hamil- 
tonian is apparent in the last term, which is proportional 
to the departure of the quasiparticle phase space distri- 
bution function n p (r, t) from the local equilibrium distri- 
bution n (ep + j): 



6n p (r, t) = n p (r, t) - n (e p+A ) 



n p (f,t) - n (e p ) - n' (e p ) 



p-A(r,t) 



n\ p - n' (e p ) 



P-A(r,t) 



(3.2) 



(2.25) Here n {e p ) = 9( 

e p F ~~ e p) i s the true equilibrium dis- 
tribution function at T = and chemical potential e PF , 
n o( £ p) = ~$( e PF ~ e p) i s ^ s derivative with respect to 
energy, and 



ni p (f,t) = n p (f,t) - n (e p ) 



(3.3) 



is the departure of the distribution function from true 
equilibrium. 

The fact that the departu re fr om local rather than true 
equilibrium appears in Eq. (3.1) is essential to guarantee 
particle conservation and gauge invariance of the theorv. 

We now make use of the well-known Landau relationta 
between bare and effective mass in a translationally in- 
variant system 



P_ 

to to 



p' 



(3.4) 



and rewrite the effective Hamiltonian in ter ms of the de- 
parture from true equilibrium [see Eq. ( |3.2|) 1: 

H qp (r,p) = e(p) + P -A^A + V f ppl n lpl (r, t). (3.5) 

TO 

P' 

Finally, we write the classical (linearized) Liouville 
equation for the evolution of the quasiparticle distribu- 
tion function under H qp . After introducing the Fourier 
representation 



n lp (r,t) = n lp (q,Lo) 



we obtain the desired equation of motion 
(q- v p - aj)n lp (q,u) 

P 



q- v p n' {e p ) 



^2f pp >n lp/ (q,Lj) 



A{q,u) 



(3.6) 

(3.7) 
= 0. 



The current response is obtained from the quasiparticle 
distribution function via the relation 



n 



= X] —ni P (q,uj) + —Ai(q,w), (3.8) 

* — ' TO TO 

P 



where it must be noted that the bare mass, rather than 
the effective mass, enters the definition of the current 
The density response is given by ni(q, u>) = n ip(9; UJ )- 
Equation ( |3.7| ) can be solved for a given value of the 
ratio x = qvp/u, with both q and u) tending to zero. 
After setting 
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ni p (q, lo) = Up(x) ■ A{q, u>), (3.9) 
we see that II^(x) obeys the equation of motion 



Hp{x) = Rp(x) [ — + ^2 fpp'Rp'( x ) 



where 



Rp{x) = n' (e p ) 



x cos(9) 
x cos(9) — 1 ' 



(3.10) 



(3.11) 



and 9 is the angle between p and q. 

In the x — > limit [see Eq. ( |l.3| )1 we expand Rp(x) and 
Ilp(x) in a power series of x as follows 



and 



Rf{x) = — n' (e p ) y^(a;cos( 



(3.12) 



(3.13) 



Inserting these expansions in Eq. (3.10) we obtain the 
recursion relation 



n' (e p )^(cos9) n 



-n' (e P )J2fpf £ n<rWr- m , (3.14) 

p' m=0 

with nl 0) = and U { } } = -n' (e p )pcos6/m. The rele- 
vant term is the one with n = 2: 



n£ 2) 



-n'o( e p) 



P 



(cosfl) 2 



£ fpp' cos n 'o( e p' ) — cos & 



(3.15) 



This enables us to calculate the current, and hence the 
response function, exactly to o rde r x 2 — {qvF /u) 2 . More 
precisely, making use of Eq. (3.8) we obtain (from now 
on, we focus on 3D) 



XL--=£— q-Uf{x)x 2 

p 

= nq 2 g| 3/5 + 4F 2 /75 + Fp/3 
mu) 2 m 2 1 + fi/3 



(3.16) 



and 



XL 



m ' — ' »" 



77! 



nq 2 g| 1 + 2F 2 /15 + F /3 
nut 2 m 2 1 + Fi/3 



, (3.17) 



where i 7 / are the usual di- 

mensionless Landau parameters!] The following partial 
results have been used to evaluate the sums over p and 
p* in three dimensions 



and 



9 I 5 2 



(3.18) 



(3.19) 



where AT(0) = m*qp /it 2 is the three-dimensional density 
of quasiparticle states at the Fermi surface. 

A direc t comparison between Eqs. (3.16, 3.17| ) and 
Eq. ( 2.19Q yields the desired expressions for f XC L(T)(®) : 



Re/ xcL (0) 



6E F §Jb + 3^2 - 



5rj 



1 + iFx 



and 



Re/^O) 



2£ F \F 2 - \ Fl 
5n 1 + \F t 



Similar results are obtained in two dimensions: 



Re/ KcL (0) 



Re/, cT (0) 



1+^1 



~n lTpi 



(3.20) 



(3.21) 



(3.22) 



(3.23) 



Substituting in Eqs. and ( |2.24| ) we obtain the 

following expressions for the elastic constants. In three 
dimensions 



K = n 



d 2 e(n) 2nE F 1 + F a 



dn 2 



3 1 + |fi 



and 



2n£ F 1 + F 2 /5 
5 l + Fi/3' 



In two dimensions 



l + F 



and 



H = 



1 + iiq " 

ngg 1 + F 2 /2 
2 l + Fi/2' 



(3.24) 



(3.25) 



(3.26) 



(3.27) 
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These are the main results of this Section. There is 
no surprise as far as the bulk modulus is concerned: 
it is given by the standard thermodynamic expression, 
where the energy density can be calculated by Quantum 
Monte Carlo. The shear modulus, on the other hand, 
has an expression involving Landau parameters, which 
are not easily calculated from Monte Carlo simulations, 
even though some progress in this direction has recently 
been reportedly 



B. Rigorous upper bounds on the elastic constants 

In this Section we derive two exact bounds on the 
elastic constants, which follow from the Kramers-Kronig 
dispersion relations between the real and the imaginary 
parts of fxcUT)^)- The origin of these relations can be 
easily seen from the formula 

fxcL(T){u) = lim ^ (jx^r 1 ^) - Xi(T)(g 5 ^)) 



(3.28) 



which directly f ollow s from the representation ( 2.17 ) and 
the definition ( p. 20 ). Both \ and x are analytic func- 
tions of frequency in the upper half plane of this vari- 
able, and both have no zeroes in this domainJlj This 
implies that their inverses are also analytic everywhere 
in th e upp er half plane. The large frequency behavior of 
Eq. ( |3.28| ) is regular because x an d X^ have the same 
form [n/m + 0(q 2 /to 2 )] in this limit. From this, one can 
also see that the q — > limit is well behaved. 

From these considerations we conclude that the 
Kramers-Kronig relations must hold in the standard form 



Re/ : 



xcL(T) 



(uj) = Kef xcL( x)(oo) 



2 f°° ,w'Im/ Icl(T) w' , 

+ -V / dJ . ' 3.29 

n Jo K) - w 

where V denotes the Cauchy principal part. The sec- 
ond law of thermodynamics (positivity of dissipation) re- 
quires that Imxipr)(?,w) < at all positive frequencies. 



But ImxfL(?,w) = for qv F < w [see Eq. fl2.18| )] 
Thus, from Eq. ( |3.28j ), we see that 



Im/xci(T)(w) < 



(3.30) 



for all positive frequencies. It then follows from Eq. ( 3.29| ) 
that 



A-cL(T)(0) < fxcL(T)(°°)- 



(3.31) 



Recall now that the right hand side of Eq. ( 3.31 ) can be 
expressed exactly, via the first moment of the current- 
current spectral function, in terms of the expectation 
values of the kinetic and potential energy in the ground 
state (a brief derivation is given in Appendix |a|) : 



/xci(oo) = — 
In 



(fce)o) 



1 + 3/JW 



(pe) 



(3.32) 



and 



fxcr(oo) = 



2/7 



U(ke)-(ke) ) + 



/?W - 1 



(pe) 



(3.33) 



where (fee), (fce)o, and (pe) are the expectation values 
of the kinetic energy, the noninteracting kinetic energy 
and the potential energy per particle respectively, and 
j3^ = 1/2, (3^ — 1/5. These quantities can be ex- 
pressed in terms of the exchange-correlation energy per 
particle e xc (n) as follows: 



and 



(fee) — (fce)o = dn 



(pe) = -dn 1+2 ' d 



l/d+l ( e xc \ 
\ n l/d) 



,2/d 



(3.34) 



(3.35) 



where d is the dimensionality and the prime denotes dif- 
ferentiation of the function in the round brackets with 
respect to n. The function e xc (n) is given in Refs. [?]|| for 
3 and 2 dimensions respectively. 

The lu — > oo limit of the longitudinal local field fac=. 
tor was first calculated in three-dimensions by Puff.E£l 
That result, is usually referred to as the "third moment 
sum rule"Ej since it is related to the third moment of 
the dynamical structure factor - the spectral function 
of the density-density response function \- Because of 
the relation \ = Q 2 Xl/u 2 (which follows from gauge- 
invariance and from the continuity equation) the third 
moment of the density-density response function coin- 
cides with the first moment of the longitudinal current- 
current response function. The straightforward exten- 
sions to two-dimensions and to the transverse case are 
outlined in Appendix [A|. 

Combining the foregoing res ults with our expressions 
for the elastic constants [Eqs. ( 3.25| - 3.24 )] we obtain the 
rigorous inequalities 



(j, < oifinEp + n 2 f xcT (co) 



K + I 2 - ~ ) n < a^nE F + n 2 f xcL (oo), (3.37) 



(3.36) 



with the coefficients a defined after Eq. (|2.20| ). The 
usefulness of these inequalities arises from the fact that 
the quantities on the right hand sides are ground-state 
properties, which can be calculated by Quantum Monte 
Carlo. Notice that the inequalities are satisfied as strict 
equalities whenever the dissipation vanishes, i.e., when 
I m /a;cL(r) = at all frequencies. In an electron liquid, 
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this happens both in the first order approximation with 
respect to the strength of the Coulomb interaction (weak 
coupling regime), and in the strong coupling limit, when, 
the electrons are expected to form a Wigner crystaLtD 
This observa tion leads u s to suggest that the right hand 
side of Eqs. ( 3.36 ) and ( 3.37 ) may provide good approx- 
imation to the elastic constants at all coupling strengths. 



IV. APPROXIMATE EXPRESSIONS FOR THE 
SHEAR MODULUS AND VISCOSITY OF AN 
ELECTRON LIQUID 

A. High-density limit 

In the regime na d B 3> 1 where as is the Bohr radius, 
the effect of the Coulomb interaction is small, and can be 
treated by first-order perturbation theory. It is straigth- 
forward to show that the f XC L(T)( ui y s are rea l an d in- 
dependent of frequency in this approximation. This is 
because, in the limit q — > 0, the imaginary part of the 
current-current response functions arises from processes 
involving at least two electron- hole pair excitations: such 
processes are not allowed in first-order perturbation the- 
ory. The vanishing of Im/ xc x,m(u>), combined with the 
dispersion relations (3.29) implies that Re f xcLCT)^) is 
independent of frequency. Hence Eq. (3.31) holds as a 



strict equality. Since (ke) = (fee)o = -^p^Ep, {pe} 



~(3/4)e 2 k F /TT for d 



3 and (pe) ■ 
2 in the first order approximation we obtain 



(4/3)e 2 fc f /vr for 



2nE F 



ne 2 k F 
IOtt 



(3 dimensions, na^ 1) and 



nEp ne 2 kp 

H - 

2 6tt 



(4.1) 



(4.2) 



naj 



^> 1). These results can also be 



(2 dimensions, 
obtained directly from Eqs. ( ft.25| , |3.27| ) of Section III A, 
making use of the first-order expression for the Landau 
parameters, fpp — —v{p — p 1 ) /2. 

Let us now turn to the calculation of the high-density 
limit of the viscosities r\ and Our starting point is 
the second-order expression for Im/i C i(T)(w), which is 
obtained from Eq. (11) of Ref. || after replacing the re- 
sponse functions Xl,t by the noninteracting ones x°l t 
and setting the "exchange correction factor" equal to 1: 



T f ( \ f ' dU ' f 



(<z) 5 



«L(T) — ImXi (q, u') + b L(T) — Imxf (q, u') 



(4.3) 



with (a L , or, b L , b T ) equal to (23/30,8/15,8/15,2/5) 
in three dimensions and to (11/16, 9/16, 1/2, 1/2) in two 
dimensions. 

The imaginary parts of the noninteracting response 
functions Xut) a ^ smau w an d finite q are directly cal- 
culated from Eq. ( |2.1g| ): 



d n 
2Ep~ 



to 
qvp 



and 



n lo 

-IT . 

m qvp 



(4.4) 



(4.5) 



for uj/vp < q < 2qp + lu/v f ; zero otherwise. The con- 
stants (7l,7t) a re given by (ir/2, 37r/4) for d = 3 and 
(1,2) ford =2. 

From the above formulas, it is easy to see that Eq. ( |4.3| ) 
gives an infinite result, due to the divergence of the 
unscreened Coulomb interaction v(q) for q — > 0. The 
result is indeed finite if the screening of the interac- 
tion is duly taken into account. In the high-density 
limit and at low frequency this is accomplished by the 
use of the Thomas-Fermi statically screened interaction 



v(q) 



,,TF 



(q) = Aire 2 /(q 2 +q 2 F ) {q 2 F = QTtne 2 /E F ) in 



three dimensions, and v(g) 



,,TF 



(q) = 2ne-/(q + q TF ) 



(qtf = 2/ob) in two di men sions. 

The first term in Eq. ([D^), which involves the product 

of two Xl^'Sj is proportional to ui 3 , and therefore does 
not c ontribute to the viscosity coefficients [see Eqs. ( [2.25j , 
2.26)] . As for the second term, we find, after some tedious 
but straightforward calculations, 



Im/ : 



xcL(T) 



R U 



„d-lr„,T.F 



v^(q)] 2 dq, (4.6) 



where /?l(t) = (3/1287t)6l(t) for d = 3 and /3l(t) = 
(l/127r 2 )6 i(T) for d = 2. 

Evaluation of the integral and substitution in 



Eq. ( 2.25 ) leads to our result for the high-density limit 
of the shear viscosity: 



,3/2 



(4.7) 



40(a B k F ) 3 / 2 60 
in three dimensions, and 
r 2 \/2 

n ~= -r-n In — = 0.053r 2 n(- In r s - 0.65) (4.8) 
6n er s 

in two dimensions. 

It is interesting to notice that, in the same limit, the 
bulk viscosity £ vanishes, both in three and two dimen- 
sions, due to the relationship &l/^t = 2(d — 
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B. Low-density limit 



In the regime na d B <C 1 the electron liquid is strongly 
correlated (via the long-range Coulomb interaction) and 
its behavior is expected to be similar to that of a clas- 
sical Wigner crystalJlZI The elastic constants of a clas- 
sical Wig-ner crystal have been calculated by various 
authors .E§[L3 Of particular interest is the case of the 
hexagonal lattice, which is expected, to be the stable 
crystal structure in two dimensions xB The elastic prop- 
erties of this lattice are formally indistinguishable from 
those of a homogeneous and isotropic body, i.e. there 
are only .-two elastic constants K and fip and they are 
given bylia /t ~ 0.24e 2 n 3 / 2 and K = -6//. (The fact 
that K < in an electron liquid should be no cause for 
alarm because this bulk modulus enters physical prop- 
erties summed to the Fourier transform of the Coulomb 
interaction v(q), which is large and positive at long wave- 
length.) In three dimensions, the Wigner crystal has 
cubic symmetry, and anisotropic elastic constants. The 
appropriate low-density limit for the strongly correlated, 
liquid, obtained by averaging over different orientationscj 
is fi ~ 0.19e 2 n 4 / 3 and K = -(10/3)/*. 

Remarkably, we find that in this case, as well as in the 
weak coupling limit, the inequality ( 3.36| ) is obeyed as 
a strict equality, namely, substituting on the right hand 
side of Eq. ( |3.36 ) the potential energy of the Wigner crys- 
tal ( l.BjuRyd in 3D, 2.2/r s Ryd in 2D), and 

neglecting the kinetic energy, which tends to zero in the 
low density limit, one obtains the correct value of /t. This 
implies that the imaginary parts of /xc£(t)(<*>)'s vanish in 
the low-density limit. 



C. Interpolation formula 

At intermediate densities no exact results for the 
fxcLlT)^)'^ are available. A mode-coupling calcula- 
tion of these quantities for both 2-dimensional and 3- 
dimensional electron gases has recently been performed 
by Nifosi et alB Their results are expected to be an im- 
provement upon previous estimates,El at least for not too 
large coupling strengths. Unfortunately, the values of /t 
obtained from this approximate theory do not conform 
to the physical expectation that /t should reduce to the 
shear modulus of the classical Wigner crystal in the limit 
of large r s . In fact, the approximate /t is found to be- 
come negative at large r s . We believe that this should be 
regarded as a failure of the approximate theory. For this 
reason, we propose an interpolation formula for \i that 
does not suffer from this problem: it reduces to the cor- 
rect limits for high and low density, and does not differ 
substantially from the estimates of Nifosi et alu in the 
range of r s where the latter are expected to be reliable. 

Our approximate formula is 




FIG. 2. Shear modulus fi/n (full curve) as a function of 
density in 3D (main figure) and 2D (inset). Dot s are from 
Ref. ^, open squares have been obtained from Eqs. ( 3.25 - 3.27 
using estimates of the Landau parameters from Refs. 
dotted curves are the asymptotic behavi ours o f Eqs. (4.1 



and from the Wigner crystal (see Section IV B), the full curve 



is the upper bound of Eq. (3.36), and the dashed curve is the 



approximate interpolation of Eq. (4.9) 



where a and b are obtained from the low-r s limit (in 
3D, a = kp/5m and b = e 2 k F /lQir; in 2D, a = k 2 F /Am 
and b = e 2 kp/6n) and c is obtained from the high-r s 
limit (c = 0.2ARyd in 3D and c = 0.22Ryd in 2D). The 
approximate /*(r s ) is plotted in Fig. together with the 
values of /t from Ref. ||. 

An analogous fit can be performed for the shear vis- 
cosity rj, with the caveat that in this case the high-r s 
behaviour is unknown and the advantage that no sign 
problem is present in the mode-coupling computation. 
The proposed formulas are 

77 ~ (60r s - 3 / 2 + cir-7 1 + car" 2 - 73 + c^rj 1 ^ 1 n (4.10) 
where c\ = 80, C2 = —40, C3 = 62 in 3 dimensions, and 



ix/n 



+ br s 1 + (c- &)- 



20 



(4.9) 



_£. l n JL= + CQ r 2 + Cl r s 2 + c 2 r s x ' 2 + c 3 
D7T er s / 



(4.11) 

where Co = 0.25, c\ — 21, C2 = 23, C3 = 13. in 2 dimen- 
sions. The resulting functions are compared in Fig. ^ 
with the values calculated in Ref. U. 
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0.06 



0.04 



.a 



0.02 - 




FIG. 3. Shear viscosity r\ in d = 2 and d = 3 (in iHnits of n) 
from the mode-coupling calculation of Nifosi et ala ( cros ses) 
compared with the analytical expressions of Eqs. (|4.7 - 4.6 ) 
(full curves). 



D. Transverse collective mode in the electron liquid 

A transverse collisionless sound mode will exist in 
the uniform electron liquid provided that the transverse 
current-current response function has a pole close to the 
real axis. Unfortunately, because such a mode must have 
a linear dispersion of the form u> — c t g the visco-elastic 
form of the response function [Eq. ( |2.15| )] - which is valid 
for u> ^ qvp - is not applicable unless Ct vp. One 
should instead use the response function in the limit 
lu I q = c t = constant, which involves all the Landau pa- 
rameters, and is presently unknown. 

The situation becomes much more favorable in the 
limit of large r s (low density). In this limit, the transverse 
sound velocity is large compared to the Fermi velocity, 
so that the visco-elastic form of t he tr ansverse response 
function can be used. From Eq. ( 2.15| ) one can immedi- 
ately deduce the existence of a pole at 



■ V 2 

*o 1 > 

Lmn 



.12) 



where cf = u/mn. This result is independent of di- 
mensionality. Note that the linewidth of the excitation 
(nq 2 jnm) vanishes in the long wavelength limit. 

From now on let us focus on the two-dimensional elec- 
tron liquid, since this is the system which, being closer to 
Wigner crystallization, provides the best chances for the 
observation of a transverse sound mode. At low density, 



making use of Eq. (4.12) and of the low-density form the 
shear modulus, wc obtain 



4- =s 0.07r s 



(4.13) 



which grows with decreasing density, and justifies the use 
of the visco-elastic form of the response function. 

Equation (4.13) can be used to give a rough estimate 
of the minimum value of r s above which the transverse 
sound mode would be observable. Requiring c t /vp > 1 
we obtain the condition r s > 14. This minimum value of 
r s is still significantly lower than the critical r s ~ 37 for 
Wigner crystallization, as estimated from Monte Carlo 
calculations B 



V. INCLUSION OF COLLISIONS 
A. The relaxation time approximation 

Up to this point, we have neglected processes that 
limit the lifetime of a quasiparticle, such as quasiparticle- 
impurity and quasiparticle-quasiparticle collisions. 

In this Section we reinstate these processes, and study 
their effect at the phenomenological level. The starting 
point of our analysis is still the kinetic equation for the 
quasiparticle distribution function, but now we include a 
collision term: 



(q ■ v p - uj)ni p (q, uj) - q- v p n' (e p ) 



V f PP >n lp , (q, ui) H A(q, u) 

— ' TO 

P' 



= — il\n 



ipJ 



(5.1) 



where I[ni p ] is the collision integral. 

Without going into the details of the collision pro- 
cess, we shall simply assume that collisions attempt to 
restore a "locally relaxed" equilibrium distribution func- 
tion n^ = nn(e n ) + n?L, with a characteristic relaxation 



ip = n (e P 
time r, namely 



/[nip] 



(5.2) 



Equation ( p.2| ) is generally referred to as the "relax- 
ation time approximation" .113 

The locally relaxed distribution function n^(r, t) is de- 
fined as the distribution that, at any given instant, would 
be in equilibrium in the presence of appropriate scalar and 
vector potentials Vn(r,t) and An(r, t), chosen so as to 
make Eq. ( |5.l| ) obey the conservation of particle number 
and (when appropriate) particle current. Let us discuss 
this construction in some detail. 

(i) Impurity scattering. In this case collisions conserve 
the local quasiparticle number (i.e., the density), but not 
the current. Therefore, we must require 
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(5.3) 



that is, the locally relaxed distribution function must 
yield the same density as the true distribution function 



(5.4) 



This is accomplished by defining as the instanta- 
neous equilibrium distribution function in the presence 
of a scalar potential Vr such that 



Vk = X _1 (0)ni(«,«) 



(5.5) 



where x(0) is the static (cj = 0) density-density response 
function in the q — > limit. Because there are no addi- 
tional constraints, the vector potential Ar remains equal 
to the real one A. 

(ii) Quasiparticle-quasiparticle scattering. In this case, 
the collisions conserve not only the local number, but also 
the local mome ntu m, i.e., the current density. Therefore, 
in addition to (5.3), we must require 



(5.6) 



that is, the locally relaxed distribution function must 
yield the same canonical current density as the true dis- 
tribution function: 

jc(q, u)=Y2 n lp (q, ui)p/m = Y n i P (<l, u)p/m. (5.7) 

p p 

On the other hand, the full locally relaxed current density 

j R = j c + nA R /m (5.8) 

must vanish, because it is the current of a system in equi- 
librium. This fixes the value of the Ar potential as 



A R (q,u) = -(m/ri)j c (q,uj). 
The value of Vr is still given by Eq. (|5.5|) . 



(5.9) 



B. Solutions of the transport equation in the 
relaxation time approximation 

Equation ([Tl]) can be solved to yield the density- 
density (x) and the transverse current-current (xt) re- 
sponse functions of the system with collisions, in terms of 
those of the system without collisions, i.e., with 1/r set 
to zero. We describe our method of solution in Appendix 

In the case of impurity scattering (no current conser- 
vation) we obtain 



to 



i/r 



X T (g» uj + i/Tx(q,u + i/T) lu + i/r x(g,0)' 



(5.10) 



where x T (<fj w ) is the density-density response function 
including collisions, and x(<7, <*> + i/ T ) is the same quan- 
tity without collisions, but calculated at the complex fre- 
quency lu + i/r. Similarly, for the transverse current- 
current response function we obtain: 



1 



lu + i/r 



1 



(5.11) 



The longitudinal current-current response function is, 
of course, obtained from the density-density response 
function via the continuity equation relation x£(g,w) = 
(lu 2 /q 2 )x T {q, w), which continues to hold in the presence 
of scattering. We note, in passing, that the above equa- 
tions, when used to calculate the electrical conductiv- 
ity, lead to the familiar Drude formula, where r is the 
electron-impurity scattering time. 

In the case of current- conserving scattering the above 
two equations are modified as follows: 



1 



1 



X T (<7>) 



lu + i/r x{q, lu + i/r) 

i/r 1 i m lu 

lu + i/r x(q, 0) t n q 2 



(5.12) 



and 



lu + i/ 



1 



c ■ ( 5 ' 13 ) 

lu XT(q, w + i/t) lut n 



Notice the additional terms on the right hand sides of 
these equations, which guarantee current conservation. 

In both cases, we define the exchange correlation ker- 
nels fxcLCr)^) m tne presence of collision, by direct gen- 
eralization of Eq. ( 3.28| ), namely 



r 



T)(^) = I™, : 



(5.14) 



where the "reference" response function XlCt)^' 1 ^) m 
the presence of collision s, is obtained from the solution 
of the kinetic equation (5.1) with all the Landau param- 
eters set equal to zero. 

We notice that our "reference function" is not the same 
as the noninteracting response function, because it con- 
tains the relaxation time which is determined, at least 
in part, by electron-electron interactions. Thus /J c is a 
mathematical construct: its purpose is to take into ac- 
count some interaction effects which admit description 
in terms of Landau parameters. Additional interaction 
effects are phenomenologically included in the relaxation 
time, and are already contained in the reference response 
function. 

With the above definitions, the collisional exchange- 
correlation kernels are found to be related to the colli- 
sionless kernels by relationships that closely parallel the 
analogous ones for the inverse response functions: 
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/J C Z,M 



UJ + %/■ 



-Jxch(u + i/r) + 



ijr d 2 e xc (n) 
uj + i/r dn 2 



and 



JxctM = fxcT(U + i/t) 



(5.15) 



(5.16) 



Note that these formulas hold for both types of scatter- 
ing. 

In practice, under the assumption that 1/r <S Ep, one 
can approximate f xcL (r){u + i/r) - /xci,(T)M- This is 
justified because the collisionless / xc 's are smooth func- 
tions of u), which vary significantly on a scale set by the 
Fermi energy (or plasmon frequency). Therefore, the 
fractional error introduced by neglecting 1/r in the ar- 
gument of fxc' s is e xpecte d to b e of order 1/EpT <C 1. 

Equations (5.15) and (5.16) provide the basis for 
an approximation to the frequency-dependence of 
the exchange-correlation kernels, which interpolates 
smoothly between the static limit and the dynamic low- 
frequency limit across a region of width 1/r in frequency. 
The form of the dependence of f xc on the inverse scat- 
tering times shows that the following additivity property 
holds: If there are two independent scattering mecha- 
nisms operating simultaneously with relaxation times T\ 
and T2 , then their combined effect is equivalent to that of 
a single scattering mechanism with an effective relaxation 
time r e ff such that 



1 



1 

n 



1 

T2 



(5.17) 



This can be prov ed st raightf orwar dly, by applying the 
transformations ( 5.1 5| ) and (5.16) twice in succession, 
the first time with relaxation time t\ and the second time 
with relaxation time . The result is the same that one 
would obtain by applying the transformation only once, 
with relaxation time r e //. 



C. Elastic constants and viscosity in the presence of 
collisions 

In order to calculate the elastic constants and viscous 
coefficien ts in the pre sence of collisions we first substitute 
in Eqs. (5.1C) and (5.11) the long- wavelength forms of 
the collisionless response functions derived in Section II. 
These are conveniently rewritten as 



nq 



1 



K 



and 



Xt(<7>) 



n 
rn 



1 + 



jji{ui) nq 2 



AM 



nq 



(5.18) 



(5.19) 



where fi(uj) is the collisionless generalized shear modu- 
lus, and where we have taken into account the fact that, 
according to our previous discussion, the bulk viscosity 
of the electron gas vanishes, and K(uj) ~ K, at low fre- 
quency. We also need the long-wavelength form of the 
static density-density response function, which is 



X(£0)~--. 



(5.20) 



Then, Eqs. ( pTlfj| ) and ( fTTl| ) yield, in the case of impui ity 
scattering, 



1 mu)(ui + i/r) uj 



X T (<?M 



nq 2 uj + i/t 

2 _2\ ftw + i/r) 



K 



-2 ( 5 ' 21 ) 



lu + i/r n 2 



and 



1 m(o; + i/r) u> £t(ui + i/r) q 2 



xTp{q,u>) nu> uj + i/r n 2 lu 2 

(5.22) 

Again, we can neglect i/r in the argument of ft with a 
relative error of order 1/EpT. 

Next, we compare the long wavelength forms of the 
response functions J5.2l| , 5.22 ) with the ones obtained 
from generalized elasticity theory in the prese nce o f col- 
lisions. In order to do this, we return to Eq. ( [2.13 ) and 
add a relaxation term —ji(q, u>)/t on the right hand side 
in the case of impurity scattering (in the case of current- 
conserving scattering no additional term is needed). The 
elastic constants must also be modified to include the ef- 
fect of collisions: we call K T and fi T the generalized elas- 
tic constants, which depend on frequency on the scale of 
1/r. Solving the modified equations of motion we ob- 
tain the density-density and transverse current-current 
response functions in the following form 



muj(uj + i/r) 



nq" 



2\ f Y(uj) 



and 



m(u) + i/T) jl T (uj) q 2 



Xt(<TM 



(5.23) 



(5.24) 



Com parin g equations (5.23) and ( |5.24 ) to (5.21) 
and ( 5.22 ) respectively we accomplish our goal of ex- 
pressing the elastic constants K T and fi T in terms of their 
collisionless counterparts K and jl: 



A r M = 



uj + i/' 



:AM 



(5.25) 



and 
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VI. CONCLUSION 




FIG. 4. Shear modulus /j, t (full curve) and viscosity rf 
(dashed curve) in the pr esence of collisions, as functions of 
lot, from Eqs. (5.27-5.28). The dots mark the lot = and 
lot = oo asymptotic limits. 



In this paper we have shown that visco-elasticity is the 
effective theory that describes the dynamical response of 
a Fermi liquid at low temperature, long wavelength and 
low (but finite) frequency. We have presented rigorous 
results and approximate expressions for the visco-elastic 
coefficients of an electron liquid in two and three dimen- 
sions. We have also shown how quasiparticle scattering 
mechanisms can be incorporated in the effective theory 
by means of a simple relaxation time approximation. An 
interesting prediction of our work is the possibility of the 
existence of a transverse sound mode in a high purity 
two-dimensional electron liquid at large r s , before crys- 
tallization occurs. 
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K t {lo) = K(lo) = K. 



(5.26) 



It is straigthforward to verify that the same results are 
also obtained in the case of current-conserving scattering. 

From Eq. fl5.26| ) we see that the bulk modulus and 
the bulk viscosity coefficient are unaffected by collisions. 
For the shear modulus and the shear viscosity we ob- 
tain, after separating the real and the imaginary parts 
of Eq. (5.25), the following equations, accurate within 
corrections of order 1 /Eft: 



(corf 



and 



\XT 



1 



V- 



1 + {lot) 2 ' 1 + (lot) 



(5.27) 



(5.28) 



In reaching the final form of these equations, we have 
used the fact that tj/t ~ 1/EpT <C [i where r\ and /i are 
the collisionless visco sity and shea r modulus. 

Equations (5.27) and ( 5.28| ) clearly exhibit the 
crossover from hydrodynamic to dynamic regime. Their 
qualitative behavior is plotted in Figure |], and one can 
observe the opposite behaviors of rj T and /i T as func- 
tions of frequency, which is consistent with the Kramers- 
Kronig dispersion relations. 



APPENDIX A: HIGH-FREQUENCY LIMITS OF 
THE EXCHANGE-CORRELATION KERNELS 

This Appendix gives a self-contained derivation of the 
high-frequency limits of the exchange-correlation kernels, 
based on the equation of motion for the current-current 
response function. 

The current-current response function Xijitfi w ) = 
(n/m)5~j,+ Rij(q,Lo) is determined by the Fourier trans- 
form of£3 



RiMt) = -<fl(t)<[7?,i(*) s J-w(0)]) 



(Al) 



where jq = J2p{p/ m ) c l-q/2 c i>+q/2 is tne canonical cur- 
rent operator. The time derivative of Rij(q,t) evaluated 
at t = gives the first frequency moment of Rij(q,Lo), 

/Oo (JjjJ d 

Ry($»W— = -sfoV.ift), j- W (0)]) 
-oo T at t=Q 

(A2) 

which is clearly real, i.e. only the imaginary part of 
Ry(<f, oo) contributes to the frequency integration. We 
now evaluate the right hand side of Eq. (|A2j) . This gives 



dt 



([?V,i(i)»J'-«(0)]) 



[[ H Jq,i],j-q,i 



(A3) 



t=0 



where H is the Hamiltonian and all operators on the 
right hand side are evaluated at t = 0. Straightforward 
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evaluation of the commutators allows us to obtain My (q) 
in terms of the momentum distribution np — (etc-) and 
the static structure factor S(q) of the electron liquid. The 
resulting expression is, up to terms of order q 2 



The collision integral (5.2) describes relaxation to a 
local equilibrium distribution n^, which was defined as 
the equilibrium solution in the presence of appropriate 
vector and scalar potentials Ar and Vr. 

In Section V we determined the value of Ar and Vr 
from the condition that the transport equation (5.1) 



TTi 



+v(p)S(p) 



d) 



q-pp q 



p>q l 



p 1 p> 

2~ 



1-dq 2 
2 



P 
d 2 



1 (p.q) : 



obey the conservation of particle number and (when ap- 
propriate) particle current. The results of those calcula- 
tions will be crucial in the following development. 

In order to facilitate the solution of the transport equa- 
tion, it is convenient to introduce a "dynamic correction 
(A4)to the quasiparticle distribution function" , defined as fol- 
lows 



where inversion symmetry has been used to eliminate 
terms of first order in q. 

From the Kramers-Kronig relations one obtains the 
high-frequency behaviour of the real part of xy(<f, a;), 

n 1 

R- e Xii(<f ; ^ -> °o) — —Sij + —Mij(q) . (A5) 

The longitudinal and transverse components of x are 
then obtained from xl = QiXijQj an d XT = UXijtj- {t 
is a unit vector perpendicular to q — q/q; a sum ove r 
repeated indices is implied). Comparison with Eq. ( 2.19| ) 
yields 



lim f xc L(T)(u) = lim -^M L(T) (q) 



oil(t)Ef 



(A6) 



In an isotropic system, the first line in Eq. ( |A4| ) is pro- 
portional to the average kinetic energy, and the re main- 
ing t o the average potential energy, leading to Eqs. ( 3.32 - 
3.33 ) of the main text. We finally r emark th at full 
isotropy is not needed to obtain Eqs. ( |3. 32 - 3. 33 ) - in- 
deed, only averages of second- and fourth-order terms 
appear in Eq. (A4). In the case of the 2D triangular lat- 
tice, such averages are identical to those of an isotropic 
fluid ((a; 2 ) = r 2 /2, (x 4 ) = 3r 4 /8, (x 2 y 2 ) = r 4 /8), and 
therefore the same results hold in the crystal and in the 
liquid. This is not the case in any of the crystals with 
cubic symmetry, such as simple cubic in 2D and fee in 
3D. 



APPENDIX B: SOLUTION OF LANDAU 
EQUATION OF MOTION WITHIN THE RTA 

This Appendix discusses the details of the computation 
of the response functions within the RTA. As in the main 
text, we denote by x T { ( li UJ ) the response function in the 
presence of a relaxation time t, and by x(<f, to + i/r) the 
collisionless response function evalutated at the complex 
frequency u> + i/r. The dependence on q, which plays 
no significant role in the computation, will henceforth be 
dropped. 



i/r 



(Bl) 



Making use of the equilibrium condition for nf p it is a 
straightforward com put ation to verify that the collisional 
transport equation (5.1) for ri\ v is equivalent to the col- 
lisionless transport equation for n® p with a modified fre- 
quency uj + i/r and modified potentials 



A n = A- 



uj + i/r 



Ar 



and 



V D = V 



i/r 
uj + %/■ 



-Vr. 



(B2) 



(B3) 



It follows that the density and current responses = 
S P n \ v i j D = S P n i P P/ m + nfA D /m are related to the 
modified fields Ajj , Vd by the collisionless response func- 
tions evaluated at frequency ui + i/r, 



V D 



(uj + i/r)q- Ai 



(B4) 



and 



3 {D) - Xl,t(" + */t)A L d T + X l(lo + i/r) 



uj + i/r 



V, 



D 



(B5) 



In writing these equations we have used the fact that 
the mixed density-current response functions Xdc and Xcd 
are purely longitudinal, and are related to the density- 
density response function by Xcd = Xdc = c&x/q 2 = 

qxL/u- 

Finally, note that n[ D \uj) = ni{ui) — -^Zppn[ R \uj), 
and that a similar relation j^(uj) — j c {^>) ~ J+i/ T jc'{ ul ) 



holds between the canonical currents and j R as well 
as between the full currents j D and j . Since j R = 
(the full current vanishes in an equilibrium state) one 
concludes that j — . We now discuss the two different 
cases mentioned in Section |y|. 
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1. Impurity scattering 



2. Current-conserving scattering 



In this case (see Section V) Ar = A\ and Vr is 
fixed from the constraint of local density conservation, 
n[ R) (q,uj) = ni(q,u). 

a. Density excitations, scalar external potential In 
this case both A and A R vanish. Density conservation 
gives [from Eq. (5.5)] 



Vn = 



Since 



AR) 



n\, we get 



x T M 
x(o) 



m = 



_uj + i/r ( D )_uj + i/i 



i.e. 



X M = '—xiu + tM 



1 



X(w + i/r)V D 



X t (uj) i/r 



x(0) uj + i/i 



(B6) 



(B7) 



(B8) 



which, after some algebra, gives Eq. (5.10). 

b. Density excitations, longitudinal vector potential 
In this case the scalar external potential V = and there 
is a purely longitudinal external vector potential A. This 
is completely equivalent to the previous case, modulo a 
gauge transformation. We carry out the computation 
only as a check. The total density fluctuation is 



q- A(q,w) 



(B9) 



and the potentials of the fictitious system are given by 
A R = A,Vr = [xEH/x(0)]<f ■ A/uj. If follows that 



3 = Xl(u + i/t) 

which gives 
uj(u> + i/r) 



u i/r g xlM 

uj + i/r (uj + i/r) 2 uj x(0) 



i/r 



A 
(BIO) 



q 2 XL (w + i/ t) q 2 xl M U) + i/r x(0) 

(BH) 

Since xl{^ + */ r ) = (w + i/T) 2 q~ 2 x( uJ + i /r) this is 
equivalent to the previous result of Eq. ( |5.1CD . 

c. Transverse excitations The external vector poten- 
tial A is purely transverse. There is no density fluctua- 
tion, V\ = Vr = Vd = 0, and Ar — A\. Thus, from the 
combination of Eqs. (B2) and (B5) we obtain 



uj + i/r 



Xt(u + i/r)A{q,uj) 



(B12) 



which is equivalent to Eq. ( 5.11 ) 



Here Vr is fixed as in Eq. (B6) , and moreover Ar is 
given by Eq. (|j|), A R = -(m/n)j c . 

d. Density excitations, scalar external potential Here 
A = and j c — j (the full current and the canonical cur- 
rent coincide). Thus Ar = —( m/ n)j = —(mu)/nq)x T Vi, 
and Vr is the same as in Eq. (|Bq). It follows that 



-Til 



1 lo + i/t 
= X(^ + i/T) 



X T H i/r . i muj T 

1 777+, 71 1 oX 

X(U) uj + i/r t nq z 



(BIS) 
V 



which is equivalent to Eq. ( 5.12| ). The same computation 
can be done using a longitudinal vector potential, with 
the same result . 

e. Transverse excitations Here V = Vr = 0; and 
Ar = —(m/n)j c = (1 — mxri^/njA, which gives 



A D = 



i/r m ' 
-Xt(u) 



A 



(B14) 



Fina lly, since Xt{ u )A 
Eq. (ftll . 



uj + i/r uj + i/r n 

Xt{u + i/T~)An we get 
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